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The paper is devoted to the relationship between the continuous Markovian description of Levy flights 

developed previously {Lubashevsky et al, Phys. Rev. E 19 (2009) 011110, 80 (2009) 031148; Eur. Phys. 

J. B 78 (2010) 207, 82 (2011) 189) and their equivalent representation in terms of discrete steps of a 

wandering particle, a certain generalization of continuous time random walks. To simplify understanding 

•/^ the key points of the technique to be created, our consideration is confined to the one-dimensional model for 

^_^ continuous random motion of a particle with inertia. Its dynamics governed by stochastic self-acceleration 

r^ is described as motion on the phase plane {x, w} comprising the position x and velocity v — dx/dt of the 

O given particle. A notion of random walks inside a certain neighborhood C of the line v — Q (the x-axis) 

S and outside it is developed. It enables us to represent a continuous trajectory of particle motion on the 

plane {x,v} as a collection of the corresponding discrete steps. Each of these steps matches one complete 
fragment of the velocity fluctuations originating and terminating at the "boundary" of C. As demonstrated, 
the characteristic length of particle spatial displacement is mainly determined by velocity fluctuations with 
large amplitude, which endows the derived random walks along the x-axis with the characteristic properties 
of Levy flights. Using the developed classification of random trajectories a certain parameter-free core 
stochastic process is constructed. Its peculiarity is that all the characteristics of Levy flights similar to the 
'T^ exponent of the Levy scaling law are no more than the parameters of the corresponding transformation from 

'^ the particle velocity v to the related variable of the core process. In this way the previously found validity 

C of the continuous Markovian model for all the regimes of Levy flights is explained. Based on the obtained 

*^ results an efficient "single-peak" approximation is constructed. In particular, it enables us to calculate the 

I I basic characteristics of Levy flights using the probabilistic properties of extreme velocity fluctuations and 

the shape of the most probable trajectory of particle motion within such extreme fluctuations. 

(N 

^ Keywords: Levy flights, nonlinear Markovian processes, random motion trajectories, extreme fluctuations, 

^^ power-law heavy tails, time scaling law, continuous time random walks 
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1. Introduction 



During the last two decades there has been a great deal of research into Levy type stochastic processes 
in various systems (for a review see, e.g., Ref. [T]). According to the accepted classification of the Levy 
type transport phenomena ^, Levy flights are Markovian random walks characterized by the divergence of 
the second moment of walker displacement x{t), i.e., ([x(t)]^) -^ oo for any time scale t. It is caused by a 
Jh power-law asymptotics of the distribution function V{x,t). For example, in the one-dimensional case this 

distribution function exhibits the asymptotic behavior V^x^t) ~ [s(i)]"/a;^+" for x 2> x{t), where x{t) is 
the characteristic length of the walker displacements during the time interval t and the exponent a belongs 
to the interval < a < 2. The time dependence of the quantity x{t) obeys the scaling law x{t) ex t^/". 
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Levy flights are met, for instance, in the motion of tracer particles in turbulent flows [5], the diffusion of 
particles in random media 3 , human travel behavior and spreading of epidemics [4] , or economic time series 
in finance [5]. 

As far as the developed techniques of modeling such stochastic processes are concerned, worthy of mention 
are, in particular, the Langevin equation with Levy noise (see, e.g., Ref. [B]) and the corresponding Fokker- 
Planck equations [7MlO|. the description of anomalous diffusion with power-law distributions of spatial 
and temporal steps [TTl [T^], Levy flights in heterogeneous media [T51 - [r7] and in external fields [IHl [H], 
constructing the Fokker-Planck equation for Levy type processes in nonhomogeneous media [20U22j . first 
passage time analysis and escaping problem for Levy flights [25^133) . 

One of the widely used approaches to coping with Levy flights, especially in complex environment, is 
the so-called continuous time random walks (CTRW) [Ml [35]. It models, in particular, a general class of 
Levy type stochastic processes described by the fractional Fokker-Planck equation (3^ . Its pivot point is the 
representation of a stochastic process at hand as a collection of random jumps (steps) {5x, 5t} of a wandering 
particle in space and time as well. In the frameworks of the coupled CTRW the particle is assumed to move 
uniformly along a straight line connecting the initial and terminal points of one step. In this case the discrete 
collection of steps is converted into a continuous trajectory and the velocity v = 6x/dt of motion within 
one step is introduced. As a result, the given stochastic process is described by the probabilistic properties, 
e.g., of the collection of random variables {v, 5t}. A more detailed description of particle motion lies beyond 
the CTRW model. Unfortunately, for Levy flights flne details of the particle motion within one step can be 
important especially in heterogeneous media or systems with boundaries because of the divergence of the 
moment /[(5x((5t)]^V Broadly speaking, it is due to a Levy particle being able to jump over a long distance 
during a short time. The fact that Levy flights can exhibit nontrivial properties on scales of one step was 
demonstrated in Ref. [301 studied the first passage time problem for Levy flights based on the leapover 
statistics. 

Previously a new approach to tackling this problem was proposed [37^0] . It is based on the following 
nonlinear stochastic differential equation with white noise £^{t) 

r^^-^v + ^r{vl + v^)m (1) 

governing random motion of a particle wandering, e.g., in the one-dimensional space M.^. Here v ~ dx/dt 
is the particle velocity, the time scale r characterizes the delay in variations of the particle velocity which 
is caused by the particle inertia, A is a dimensionless friction coefficient. The parameter Va quantifies the 
relative contribution of the additive component ^a (t) of the Langevin force with respect to the multiplicative 
one w^m(i) which are combined within one term 

VaUt) + <™(t) ^ ^{vl^V^)m . (2) 

Here we have not specified the type of the stochastic differential equation (IT]) because in the given case all 
the types are interrelated via the renormalization of the friction coefficient A. It should be noted that models 
similar to Eq. (fTl) within replacement pi can be classified as the generalized Cauchy stochastic process [41] 
and has been employed to study stochastic behavior of various nonequilibrium systems, in particular, lasers 
HI], on-off intermittency [43:, economic activity 03], passive scalar field advected by fiuid [45 , etc. 

Model (IT]) generates continuous Markovian trajectories obeying the Levy statistics on time scales t^r 
Using a special singular perturbation technique [38 it was rigorously proved for the superdiffusive 
regime matching 1 < a < 2 [37l [38] and also verified numerically for the quasiballistic (a = 1) and 
superballistic (0 < a < 1) regimes [39]. Moreover, the main expressions obtained for the distribution 
function V{x,t) and the scaling law x{t) within the interval 1 < a < 2 turn out to hold also for the whole 
region < a < 2 [S^- After its generalization [40^ model (IT]) generates truncated Levy flights as well. 

The given approach can be regarded as a continuous Markovian realization of Levy flights. Indeed, 
it is possible to choose the system parameters in such a way that, on one hand, the "microscopic" time 
scale T be equal to an arbitrary small value given beforehand and, on the other hand, the system behavior 
remain unchanged on scales t ^ t [371 [SB]- The goal of this paper is to elucidate the fundamental features 



of this approach, to explain the found vahdity of model (fTl) for describing Levy flights of all the regimes, 
and to construct a certain generalization of the continuous time random walks that admits a rather fine 
representation of the particle motion within one step. The latter feature will demonstrate us a way to 
overcome the basic drawback of the CTRW model, it appeals to different mechanisms in modeling the 
particle motion within individual steps and the relative orientation of neighboring steps. 

2. Continuous Markov model of Levy flights 

2.1. Model 

Following [371 [321 HO] let us consider random walks {x{t)} of an inertial particle wandering in the one- 
dimensional space Ka; whose velocity v — dx/dt is governed by the following stochastic differential equation 
written in the dimensionless form 

^ = -avk{v) + V2g{v) o ^{t) . (3) 

Here the constant a is a system parameter meeting the inequality 

< a < 2 , (4) 

the positive function k{v) > allows for nonlinear friction effects, ^{t) is the white Gaussian noise with the 
correlation function 

{mm)^s{t~t'), (5) 

and whose intensity g{v) > depending on the magnitude the particle velocity describes the cumulative 
effect of the additive and multiplicative components of the Langevin random force as noted in Introduction. 
The coefficient ^/2 has been introduced for the sake of convenience. The units of the spatial and temporal 
scales have been chosen in such a manner that the equalities 

fc(0) = 1 and 5(0) = 1 (6) 

hold. Equation ([3| is written in the Stratonovich form, which is indicated with the multiplication symbol 
o in the product of the white noise ^(t) and its intensity g(v). To avoid possible misunderstanding we note 
that this equation has been written in the Hanggi-Klimontovich form in Refs. [33 [311 [1^. In the kinetic 
theory of gases quantities relative to k{v) and g{v) and describing regular and random effects caused by the 
scattering of atoms or molecules are usually called "kinetic" coefficients (see, e.g., Ref. [46 ). For this reason 
the functions k(v) and g{v) together with other c^uantities to be derived from them will be also referred to 
as the kinetic coefficients. 

In the present paper the main attention is focused on the special case when the kinetic coefficients k(v) 
and g{v) are of the form 

fco(v) = 1 and go{v) = V^ + v^ (7) 

and the random walks {x(t)} generated by model ^ can be classified as Levy flights for large time scales, 
i.e., for t ^ 1 in the chosen units [37H32]- However, where appropriate, the general form of these kinetic 
coefficients will be used to demonstrate a certain universality of the results to be obtained and the feasibility 
of their generalization, for example, to the truncated Levy flights [40]. Nevertheless, in further mathematical 
manipulations leading to particular results the following assumption about the behavior of the kinetic 
coefficients 

? ' ft* I 7? 1 

k{v) w 1 , g{v) WW for 1 ^ w < Vc , and > B for v > Vc (8) 

will be accepted beforehand. Here 5 > is some positive constant and Wc 3> 1 is a certain critical 
velocity characterizing the region where the generated random walks deviate from Levy flights in properties. 
Naturally, case ((tI) obeys these conditions in the limit Vc -^ oo. 



To make it easier to compare the results to obtained in the following sections with the characteristic 
properties of the velocity fluctuations governed by Eq. (Is]), here let us find the stationary distribution 
function P^^{v) of the particle velocity v. The Fokker-Planck equation describing the dynamics of the 
velocity distribution P{v,t) and matching Eq. ([3| is written as (see, e.g., Ref. [17]) 



dP d ( . ,d[g{v)P] , , ,„ 



Its stationary solution P^'(u) meets the equality 



whence it follows that 



giv) 



d[g{v)P' 
dv 



P'\v) 



5(^') 



exp 



avk{v)P^ 



uk{u) 
-a I ,,, s du 

9 m 



(9) 



(10) 



where the constant 



C'y 



dv 
9{v) 



exp 



ukiu) , 
-a I ^, . du 
9 m 



is specified by the normalization of the distribution function P{v,t) to unity. 



(11) 



2.2. Additive noise representation 

To elucidate the general mechanism responsible for anomalous properties of the stochastic process at 
hand let us pass from the velocity f to a new variable rj — ri{v) introduced via the equation 



1 



drj 

dv g{v) 



subject to the condition r/(0) — . 



(12) 



Equation (Is]) is of the Stratonovich form, thus, we may use the standard rules of change of variables in 
operating with it. Therefore relation (12) between the variables v and r] enables us to reduce Eq. pi to the 
following one 

J ^ -acj^irj) + V2m , (13) 

where the function (j){r]) is specified by the expression 

vk{v) 



Hv) 



9{v) 



(14) 



D— ii(r^) 



The function (piv) together with the parameter a determine the rate of the particle regular drift in the 
?7-space, M,,. In particular, in case ([7]) we have 



V = sinh(77) and </'o('7) — tanh(77) . 



(15) 



It should be noted that the multiplication symbol o has been omitted in Eq. (13) because the noise ^(i) 
enters it additively and, thus, all the types of this stochastic differential equation have the same form (for 
details see, e.g., Ref. [47]). Let us also introduce into consideration the potential 



t)(r)) 



Hv) 



0(C)dC 



uk{u) 



du 



(16) 



which will be used below; in case ([7]) it can be written as 

$o(?7) = In [cosh(77)] 



(17) 



Equation (131 together with expression (14) enable us to find the time pattern {ri{t)} for a given realization 

of the noise lew }t;S- 

Leaving expression (12 1 aside for a moment and keeping in mind solely the governing equations ([3]) and 
(13 1 we may state that the constructed stochastic process {r]{t)} and the random variations {v{t)} of the 
particle velocity are related to each other via the white noise ^(i) and are of the same level of generality. 
Only relation ( 12 ) treated as the definition of the function rj = r](v) causes us to regard the time variations 
{77(^)1 as a stochastic process derived from {v{t)}. However, expression (12) may be read in the "opposite 
way" as the definition of the function v = v{rj). In this case the process {r]{t)} plays the role of the noise 
source endowing the particle motion with stochastic properties and the particle velocity v = vijf) becomes 
a derivative random variable. The particle displacement x{t) in the space M^; during the time interval t 



x{t) = / V [r^{t')] dt' 



(18) 



is also a derivative stochastic process of {'ri{t)}. Integral (18) can be treated in the Riemann sense for a 



given pattern {?y(i)} because the correlation function {ri{t)r](t')) of the stochastic process ri{t) is smooth, in 
particular, has no singularity at i = t'. It should be noted that the situation would be much more complex if 
we dealt with nonlinear integrals of white noise because the correlation function of white noise proportional 
to S{t — t') comes up at the boundary of integration regions [47] . 

As in the previous section let us find the stationary distribution function p^^{ri) of the variable ry which 
will be used below. Going in the same, way we appeal to the Fokker-Planck equation 



dp d 
dt di] 



dp d^iri) 

IT +"^^ — P 
or] arj 



(19) 



matching Eq. (13) and governing the evolution of the distribution function p{ri,t). In the limit i — >■ cx) its 
solution gives us the desired result 



f'iv) = 



2 / £-"*("') drj' 



,-Q*(77) 



(20) 



It should be noted that Exps. (10) and (20), as it must, coincide with each other within the cofactor drj/dv 
(see Exp. (12|) coming from the transformation of the elementary volume in the transition ]R„ — >■ M^. 



2.3. Trajectory classification 

As shown in Refs. [371 I38j . it is extreme fluctuations in the particle velocity, i.e., large amplitude peaks 
in the time pattern {v(t)} that are responsible for the random walks at hand exhibiting the properties of 
Levy flights. Figure [I] illustrates the multiscale structure of the pattern {v{t)}. It turns out that for any 
time interval i ^ 1 only a few peaks with the largest amplitudes inside it contribute mainly to the particle 
displacement j37j . The amplitudes of such velocity fluctuations are distributed according to the power-law 
and their characteristic magnitude grows with the observation time interval t also according to the power- 
low. The former feature gives rise to the Levy distribution of the particle displacement, the latter one is 
responsible for the Levy time scaling [38 . 

These peculiarities make it attractive to employ the following classiflcation of particle trajectories. The 
general idea is to partition any random trajectory {x(i)} in such a manner that its fragments {x{t)Yi 
or, more strictly, their "projections" {w(t)}*^^ onto the velocity space My can be treated as random walks 
inside or outside a certain neighborhood C of the origin w = 0. In this case it would be possible to 




Figure 1: The typical form of the time pattern {v{t)} generated by model (jsl in the frameworks of ansatz ^. The two frames 
depict the same time pattern for different time scales. Based on the results presented in Ref. |37| . the value a = 1.6 was used. 



Diagrams of possible trajectories: 
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Figure 2: The main diagrams and their basic fragments of the random walk classification enabling us to represent a particle 
trajectory in the space R,, as a countable sequence of alternate fragments matching the particle motion "inside" and "outside" 
a neighborhood C of the origin 77 = 0. For the sake of simplicity the figure shows only trajectories starting from the point rjy^ 
and getting a certain point tj at time t. 



regard the fragments of random walks outside the region H, as individual peaks of the time pattern {u(i)}. 
Naturally, instead of the space M^ the space M,, of the variable 77 can be equivalently used, so in what follows 
the corresponding neighborhood in the space M^ will be designated with the same symbol £ to simplify 
notations. Besides, to operate with the desired partition efficiently it has to be countable. 

A direct implementation of this idea, however, faces a serious obstacle. The matter is that a random 
trajectory is not a smooth curve. Therefore, although for a wandering particle it is possible to calculate the 
probability of getting the boundary of £ for the first time, the question about crossing this boundary for the 
second time is meaningless. The particle will cross the boundary immediately after the first one and ordering 
such intersection points in the countable manner seems to be just impossible. The trajectory classification 
to be constructed below enables us to overcome this problem and, thus, to implement the desired partition. 
Previously the key elements of this classification were used in describing grain boundary diffusion [48] and 
subdiffusion along a comb structure [15] . 

Figure [2] demonstrates the required classification of the particle random walks {?7(t)} in the diagram 
form. Let us describe it step by step. 

(js) First of all, just to simplify subsequent mathematical manipulations we confine our consideration to 
the upper half-space M+ = {?; > 0} assuming its boundary, 77 — 0, to be reflecting. The symmetry of the 
kinetic coefficient 0(77), namely, </>(— 7/) = 't'ij]), allows it. Then two values r\i and t]u such that 

< 77, < 7/„ (21) 

have to be introduced. The choice of their specific magnitudes will be reasoned below. Here, leaping ahead, 
we note that in the frameworks of ansatz ([7|) the value 77, should not be a too large number enabling us, 
nevertheless, to approximate the kinetic coefficient (j^oiv) = tanh(77) by the sign-function for I77I > 77/, 

1 if 7/ > , 

0OW « sign(77) = <( if7/ = 0, (22) 

1 if7/<0. 



For example, for 77, = 2.0 the estimates of tanh(77/) w 0.96 and vi = sinh(7;;) w 3.6 hold by virtue of (15). 
The value 77^ just should not exceed rji substantially. So in some sense 

1 < 77, and 77„ - 77, < 1 (23) 

is the optimal choice. Certainly, the final results describing the particle displacement in the real space Mj. 
do not depend on the two values. 

(js) Again to simplify mathematical constructions without loss of generality let us consider trajectories 
{'7(^')}t'=o originating from the point 77^, i.e., set r]{t')\t'=o — rju- Their terminal point 77 — ri{t')\t'=t can 
take an arbitrary value. 

(ss) Now we are able to specify the desired partition. Each trajectory {'rj{t')Yt'Z!o i^ represented as a 
certain collection of fragments shown in Fig. [2] which is written symbolically as 

G^") (r;) = ¥°f ® Fi"2 ® V°f ® F!j"4 ® . . . ® G°"*/'"(77) . (24) 

If there is only one fragment in the collection Gj" (77), namely, the last fragment, then the number n of the 
internal fragments is set equal to zero, n = 0. The meaning of these fragments is as follows. 

o The first fragment Fg"* represents the motion of the particle outside the interval £1 = [0,77,) until it 
gets the point 77, for the first time at a time moment ti > 0. Further such particle motion will be 
referred to as random walks outside the neighborhood C of the origin 77 = 0. 



o 



The next fragment F"2 matches the particle wandering inside the interval £„ = [0, 77^) until it gets the 
point 77„ for the first time at a certain moment ^2 > ^i. Particle motion of this type will be referred to 
as random walks inside the neighborhood C. 



o 



The following two fragments F23' and F34 are similar to the first and second ones, Fg"* and F'/^, 
respectively. The time moments corresponding to their terminal points are designated as ^3 and t^. A 
sequence of such fragments alternating one another makes up the internal part of the given trajectory. 



Sequence (24) ends with a final fragment which can be of two configurations, G° (77) and G'"(77). 



- The former one, G°"'(?7), is related to the particle motion starting from the point rju at time i„ 
and getting the point rj at time t without touching the boundary point rj = rji oi the interval 
Ci. This configuration exists only for rj > -qi and is characterized by an even number of the 
intermediate points, n — 2k, of the trajectory partition. 

- The latter one, 6^(77), is similar to the former configuration within the exchange of the start and 
boundary points; now rji is the start point whereas ?7„ is the boundary point of the interval £„. 
For the configuration G™(r]) to exist the terminal point 77 has to meet the inequality "q < rj^. The 
corresponding number of the partition points is odd, n — 2k + 1. 

It should be noted that when the terminal point 77 of the given trajectory belongs to the interval 
rji < q < rju both the configurations exist. 

o In addition, for rj > rji there is a special configuration G^ (77) of the trajectories that start from the 
point 77„ at the initial time t = and get the point 77 at time t without touching the boundary point 



of £1, i.e., 77 = 77;. It represents sequence (24) with no internal fragments, 77, = 0. This configuration is 
actually equivalent to G°"'(7;) with i„ — 0. 



(js) The introduced fragments of particle trajectories may be characterized by other classification pa- 

..out/in 



rameters, denoted symbolically as ©°" /™, in addition to the initial and terminal time moments of their 



realization, ti and t^+i. Figure p^ shows such parameters for the terminal fragments of Gj (77) too. The 

collection {0°"+i"} i^ determined by the specific details we want to know about the time pattern {ri{t')}l,Zo- 
In the present paper all the random walks in M„ will be separated into different groups considered individu- 
ally according to the largest amplitude 9 attained by the variable t] during the corresponding time interval 
[0, i]. As will be clear below, this classification can implemented imposing the addition requirement on each 
fragment of random walks outside the region C that bounds variations of the variable t] inside it. Namely, 
we assume the variable 77 not to exceed 9, 

■qi < T]{t) < 9 for ti < t < ti+i and Vi . (25) 

In this case each set ©H+i = {^} contains only one parameter taking the same value for all of them. For 
the terminal fragment G°"'(77) a similar condition holds. For the random walks inside the region C there 
are no additional classification parameters, so the sets ©J^+i are empty. ■ 

The constructed partition enables us to develop a more sophisticated description of the pattern {77 (t)}, 
in particular, consider its configurations where the amplitudes of its m largest peaks take given values 
9i,92, ■ ■ ■ ,9m- For sure, such analysis will allow us to penetrate much deeper into the properties of Levy 
fiights, which, however, is worthy of individual investigation. 

The remaining part of the paper will be devoted to the statistical properties of random walks described 
in terms of the constructed partition. Before this, however, let us discuss the relationship between the given 
classification of the particle motion and the well-known model of continuous time random walks (CTRW). 

2.4- Constructed partition as a generalized CTRW model 

The CTRW model imitates random motion of particles by assigning to each jump of a wandering particle 
a jump length x and a waiting time t elapsing between two successive jumps which are distributed with a 
certain joint probability density 7/^(0;, t). When it is possible to write this probability density as the product of 
the individual probability densities of jump length, ipx{x), and waiting time, tptit), i.e., ipi^i t) = '4'x{x)'ipt{t), 
the two quantities can be regarded as independent random variables and the model is called decoupled 



CTRW. Another widely met version of this model, coupled CTRW, writes the probability density ip{x, t) as 
the product of two functions 'ip{x^t) = '4>j:{x)'4>v{x/t) or 'ip{x,t) = 'ipt{t)^v{x/t). In the given case the jump 
parameters x and t are no longer independent random variables; the independent variables are the jump 
length X (or the waiting time t) and the mean particle velocity v ~ x/t. The continuous implementation of 
the coupled CTRW is based on the assumption that within one jump the particle moves along the straight 
line connecting its initial and terminal points with the fixed velocity v = x/t. 

The constructed partition can be regarded as a certain generalization of CTRW that allows one to 
consider the detailed structure of elementary steps. Indeed, a pair of succeeding fragments of random walks 
inside, ^T.i+n ^ii'i outside, F°"\ j_,_2, the region C form an elementary step of the equivalent discrete random 
walks whose jump length and waiting time are random variables. If the random walks inside the region C do 
not contribute substantially to the particle displacement x we may speak about a stochastic process similar 
to coupled CTRW. The main difference between the two processes is due to the fact that the particle under 
consideration does not move uniformly within the fragment F°"\ ^^3 • So the knowledge of the corresponding 
mean velocity is not enough to describe the real particle motion. The random walks inside the region C 
can be essential for the particle motion in space when, for example, a wandering particle spends the main 
time inside the region C. In this case the corresponding stochastic process may be categorized as coupled- 
decoupled CTRW. In the case under consideration, as will be demonstrated below, the random walks inside 
the region £ are not responsible for the Levy type behavior of the particle motion. However, the constructed 
partition actually does not require the normal behavior of random walks near the origin w = 0, so it could 
be also used for modeling stochastic processes with kinetic coefficients exhibiting singularities at t; = 0. 
Under such conditions the random walks inside the region C are also able to cause anomalous properties of 
wandering particles, giving rise to the power-law distribution of the waiting time. 

2. 5. Green function 



Continuing the constructions of Sec. 2.3 the present section analyzes the statistical properties of the 
random trajectories {??(t')}i'=o meeting condition (25) that are imposed upon all the fragments of the 
random walks outside the region C It will be done based on the calculation of the Green function Q{ri,t), 
i.e., the probability density of finding the particle at the point rj at time t provided initially, i = 0, it is 
located at the point rju and does not cross the boundary rj — within the time interval (0, t). The developed 
classification of random walks illustrated by the diagrams in Fig. [2] enables us to write 

Giv, t) = eH(ry - m)Q^ iv, t) + ©h(^« - v)G^ iv, t) , (26) 

where we have introduced the functions 

00 „ 
GHv,t) = Y. dt'g°^\7^,t- muit') , (27a) 

t t' 
00 „ „ 

G\v.t) = E / ^^' / dt"g''\il.t- t')F°-\t' - t")^k{t") , (27b) 




eH(r?)= ' . ' ' (28) 



and the Heavisidc step function 



to combine the cases < rj < rji, rji < rj < rju, and rju < rj < 9 in one formula. Here 5°"*(?7,t) is the 
probability density of finding the particle at the point rj at time t provided it starts from the point 77^ 
remains inside the region rn < -q < 6 within the time interval Q <t' <t. The function Q™{ri,t) is actually 
the same probability density within the replacement of the start point rju — > rji and now the localization 
region is < 77 < rj^. The function ^fc(i) specifies the probabilistic weight of the pattern P(i|fc) shown in 
Fig. Is] which comprises k fragments of random walks outside the region C and k fragments of random walks 



Diagram of the peak pattern P(t | k) 
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Figure 3: Diagram illustrating the structure of the pealt pattern P(i | k) for fc > 0. 

inside it. The former fragments will be referred also as to peaks of the pattern P(t|fc). The probabilistic 
weight of the pattern P(t|fc) is given by the following formula for /c > 



*PfeW = 



0<ti<t2---<t2k-l<t 

with to = and for k — 0, by definition 

%W - Sit) . 



k-l 

dtidt2 ■ ■ ■ dt2k-l _[ _[ -7^'"(^2i+2 — ^2i+l)-^°" (^21+1 " t2i) 
1=0 



(29a) 



(29b) 



In Exps. (27b) and (29a I the functions J-"°"*(t) and J-"™(t) with the corresponding values of the argument 
t determine the probabilistic weights of the fragments F°"^]^ and ¥f\_^_^, respectively. The meaning of the 
former function is the probability density of reaching the point iji for the first time at the moment t provided 
the particle is located initially at the point 77^. The meaning of the latter function is the same within the 
replacement ??; O 77^. To simplify the notations, the quantities 77^, 77;, and 6 have been omitted from the 
argument list of the corresponding functions noted above. 



Integrals (27) and (29a) are of the convolution type, therefore, it is convenient to pass from the time 



variable t to the corresponding variable s using the Laplace transform 



W{s) := W^(s) = / dte-''W{t) 
Jo 



(30) 



where Wit) is one of the functions entering Exps. (27) and (29). In these terms the given integrals are 
reduced to algebraic expressions, namely. 



and, thus. 



-piis)=[F-^^s)F°-'is)]' 



g'^('m) = E [^"(^)^™*(^)]'g°"'(^' 



G°"*(77,s) 



CXD 



G"('7,s) = 5] [F-(s)F°"*(s)]V™'(s)G-(77,s) = - 



1-F'n(s)i^°"t(s) ' 

F°"'(s)G'"(7;,s) 



fc=0 



-F'"(s)F°'-"(s) 



(31) 

(32a) 
(32b) 



because |F'"(s)F°"*(s)| < 1. 
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As demonstrated in Appendix A the functions F'"(s), F°"'(s) sa 1 for s ^ 1 provided, in addition, the 



inequahty 9^1 holds. Therefore dealing with time scales i ^ 1 only the dependence of the denominator on 



the argument s should be taken into account in evaluating Exps. (32 1. The other terms may be calculated 



within the limit s — > 0. In this way, using formula (A. 26b) derived in Appendix A. 2 the Green function 
components G'^'^ {q, s) determined by Exps. (32) can be represented as 



G^^^V.s) 



dp'^T]') 



drj' 



aG° 



'{ri.s) 



[e<^Hv.) - e"*(')0] /^°° d7y'e-"*("') 



(33) 



its Laplace transform G(?7, s) specified by Exp. (26 1 is written as 



in the leading order. Then appealing to Exp. (A.6b) the desired Green function or, speaking more strictly. 



G(r7,s) 



dp^^ (rj' 



dri' 



p-'iv)- 



(34) 



It should be reminded that here p^^ (rj) is the stationary distribution function of the random variable rj after 
merging the half-spaces {rj > 0} and {77 < 0}. 



Expression ( 34 ) for the Green function together with the constructed probabilistic weight ^k (t) of the 
peak pattern P{t\k), Exps. (29) and ([sl]), are the main results of the given section. They enable us to 



consider the probabilistic properties of individual peaks in the pattern P(i|/c) instead of random trajectories 
as whole entities, which is the subject of the next section. 



3. Peak pattern statistics 

As follows from the governing equation ^ fluctuations of the particle velocity v are correlated on time 
scales about unity, t ^ 1. Therefore, in the analysis of Levy flights based on model ([3]) within the limit 
t 3> 1 the information about the particle velocity 77 (or v in the initial units) is redundant. So to rid our 
constructions of such details let us integrate Exp. (34) over all the possible values of the velocity rj, i.e. over 



the region < r] < 9. Since in this context the influence of the upper boundary ^ 1 is ignorable, the 
integration region may be extended over the whole interval < 77 < 00. Then denoting the result of this 



action on the left-hand side of Exp. (34) by 



©l(s)= / Gir^,s)drj 



(35) 



we get 



e5i(s) 



dp^^{ri) 



drj 



ri=e 



r£[F'"(5)F°"*(5)]'. 



(36) 



k=0 



In deriving the second equality of Exp. ( |36[ ) formula (A. 26b) for the probabilistic weight of the composed 
unit F°"*^j^ (E) ¥\"^-^ ^^2 of the peak pattern P(< | k) has been used. Namely, first, it has enabled us to write 



F'"(s)F°"'^(s) = l-r 



where the introduced time scale 



dp^^ (rj) 



dr] 



Ti=e 



(37) 



ga$(rj„) _ ga*()7;) 



drj e 



/ -Q$(j)') 



(38) 
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can be interpreted as the mean duration of this unit. FinaUy, the equality 






fc=0 



F'"(s)F°"t(s) 



has led us to Exp. (36) 



Since the hmit of large time scales, t ^ 1, is under consideration, the inequalities s <C 1 and 6^1 are 



assumed to hold beforehand. In this case, by virtue of Exp. (37) 

l-F'"(s)F°"'(s)<l 



and in sum ( 36 ) the terms with A: ^ 1 contribute mainly to its value. Thereby we may confine our 



consideration to peak patterns ¥{t \ k) that are composed of many elementary units. It allows us to convert 
the discrete sum X]fc=o(- ■ ■) ^^^^ ^ continuous integral as follows 



oo 
oo ^ 

E(---)^ /dfc(...) 

fe=0 i 



and make use of the approximation 



[i^''^(s)i^°"'(s)] =exp<^-fcr 



dp^^ (rj) 



drj 



n=e 



In these terms formula ( 36 ) reads 



&l{s) = t dkcxp ■{ —kr 





dp'\v) 



drj 



ri=e 



oo oo ^ 

T dt dke^'*6{t - kr) exp I kr 

n n ^ 



dp^^ {rj) 







drj 



(39) 



(40) 



(41) 



ri=e 



The last expression enables us to represent the original function ©(i) of the Laplace transform ©l(s) as 

dp^\r,) 



(3{t) ^T / dkS{t- kr) exp { kr 



drj 



Uexp^^ 



(42) 



ri=S , 



In other words, within the given limit it is possible to consider that only the pattern P(t | k) with k = t/r 
(or more strictly k = [t/r]) contributes to the function ©(i). 

The latter statement is the pivot point in our constructions that admits another interpretation of the 



us rewrite Exp. (42) in the form 



function (S{t) and its efficient application in the theory of Levy flights. Appealing to approximation (40) let 

t 



<8{t) = [r"(0)F°"*(0)] 



fct 



with 



h 



(43) 



In what follows we will not distinguish between the ratio t/r and the derived integer fcf = [t/r] because for 
fc 3> 1 their difference is of minor importance and will keep in mind this integer where appropriate in using 



the ratio t/r. Therefore formula (43) can be read as 



oo 

e{t) = f^k,{t')dt' for t » T ~ 1 



(44a) 



or, by virtue of ( 29 ) , 










(44b) 



0<tl<t2---<t2fc<OO 
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Expressions ( 44 ) are one of the main technical resuhs obtained in the present work. In particular, they allow 



us to state the following. 

Proposition 1. The probability &{t) of finding the particle whose "velocity" rj is localized inside the region 
|?7| < 6 during the time interval (0,i) for 6*,^ ^ 1 is equal to the probability *Pfcj(i') of the pattern V{t'\kt) 
with kt = t/r peaks after the averaging over its duration t' . 

It should be noted that the given statement has been directly justified assuming the initial particle velocity 
to be equal to r/g = t]u- However, as can be demonstrated, it holds for any rjQ <C 0. Also it is worthwhile 
to underline the fact that within this description the physical time t is replaced by the number kt of peaks 
forming the pattern P(i'|fcj), whereas its actual duration t' does not matter. It becomes possible due to a 



certain self-averaging effect. In particular, according to Exp. (371, the duration of the units {F°"lj_i(8)F" ^ ^^j} 
should be distributed with the exponential law, at least, on scales t ^ t. As a result the duration t' of the 
pattern V{t'\kt) changes near its mean value fc^r = t with relatively small amplitude. Besides, to simplify 
further explanations, when appropriate the probability &{t) will be also referred to as the probability of the 
pattern P(fct), where the argument t', its duration, is omitted to underline that the required averaging over 
t' has been performed. 



Expression ( 44b ) admits another interpretation of the developed random walk classification. 



Proposition 2. The probability 25 (t) is equal to the probability of the particle returning to the initial point 
V — Vu for kt — t/r times sometime after. The return events are understood in the sense determined by the 
sequence of alternate jumps between the two boundaries {?7;,?7m} of the layer C. 

The following proposition will play a significant role in further constructions. 

Proposition 3. The probability ©(t) of the peak pattern P(fct) is specified by the expression 



e(t) = exp <^ t ^PWZ 
drj 



(45) 



and the probability of its one unit, i.e., the probability of one return to the initial point rj = rju sometime 
after is evaluated by the expression 

0(t) = l + T ^ ^'^ 



drj 



(46) 



Formula (45) stems directly from Exp. (41), whereas formula (46 1 is a consequence of Exps. (43) and (37) 



where in the latter one the parameter s is set equal to zero, s = 0. 

4. Core stochastic process 

As noted above spatial displacement of the wandering particle is manly caused by extreme fluctuations 
in its velocity. Therefore the direct contribution of the random walks inside the layer C to the particle 
displacement is ignorable. We have to take into account only the fact that during a certain time the 
wandering particle or, more strictly, its velocity rj (or v in the original units) is located inside this layer. 
Thus, if two given stochastic processes differ from each other in their properties only inside the layer C 
but are characterized by the same probabilistic weight jF™{t) of the fragments {Fi"-|_i}, then they can be 
regarded as equivalent. In the case under consideration, i.e., for i 3> 1 it is sufficient to impose the latter 
requirement on the Laplace transform F^'^{s) of this weight within the linear approximation in s <C 1. 

The purpose of the present section is to construct a certain stochastic process that, on one hand, is 
equivalent in this sense to the process governed by model ([3]). On the other hand, its description should 
admit an efficient scaling of the corresponding governing equation, which will enable us to single out the 
basic mechanism responsible for the found basic features of particle motion. The other characteristics of 
particle motion can be explained by appealing to the corresponding coefficients of this scaling. The desired 
equivalent process will be called the core stochastic process for these reasons. 
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Figure 4: Illustration of constructing an effective stochastic process {»?'(<)} with the V-type potential that is equivalent to the 
initial process {»?(*)} from the standpoint of the particle displacement in the space R^. 



4.I. Construction 

The key points in constructing the core stochastic process are illustrated in Fig. |4] The construction is 
based on the transformation 



ri = ri' + sign(77')Aa 
and the replacement of the potential $o {v) by its linear extrapolation to the region "q <1 



^oiv) 



\ri 



A,- 



ln2. 



(47) 



(48) 



The constant Aq, has to be chosen such that in both the cases the function F™(s) take the same form 
within the linear approximation in s. Appealing to Appendix A. 2.1 namely, formula (A. 11) we see that 
this requirement is reduced practically to the equality 



-a*o(i;) 



drj 



-Q!(r;-ln2) 



drj 



(49) 



because the potential <I>o (?]) differs considerably from its linear interpolation ( 48 ) only in the region rj <\ 



The direct calculation of the integrals entering equality (491 for function (17 1 yields the desired value 



A„ 



In 



4r(a) 



aV^{al2) 



(50) 



where r(. . .) is the gamma function. Figure Is] shows the value Aq, as a function of the parameter a € (0, 2). 
It should be noted that previously the value of Aq, has been implicitly assumed to be less than rji , which is 
justified by Fig. [5] Indeed, on one side, the boundary r^i is initially assumed to belong to the region wherein 
the function $o('7) admits the linear approximation, i.e., the estimate ^7; ^ 1 is assumed to hold beforehand. 
On the other side, the maximal value of A^ is about 0.35 according to Fig. [5] 

In what follows we will confine our consideration to the special case (IT]) matching the ideal Levy flights. 
It enables us to convert from the stochastic process {ri{t)} to the effective stochastic process {r]'{t)} using 
transformation (47 1 and to reduce Eq. (13) to the stochastic equation 



di 



-asign{T]') + V2£,{t) 



(51) 
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Figure 5: The magnitude of the parameter A^ vs the possible values of the parameter a. 



governing the analyzed Levy flights in the equivalent way. Then, via the next transformation 77' — > u specified 
by the expressions 

?7' = -u, i=^t (52) 



involving also the transformation of the time scales t — > t equation (51 1 is rewritten as follows 



di 



- sign(u) + ^C(t) . 



It is a parameter-free stochastic differential equation with additive white noise such that 

(e(t)) = o, (c(iK(t')> = 5(t-t'). 



(53) 



(54) 



In other words, we have constructed the stochastic process {u(t)} governed by Eq. (53 1 that can be treated 
as the desired core process. It is of the same form for all the types of one-dimensional Levy flights, at 
least, Levy flights described by models similar to Eq. (Is]) inside the region where the cut-off effects are not 
significant. The basic parameters characterizing the generated Levy flights such as the exponent of the Levy 
scaling law depend on the system parameters, in particular, the coefficient a via the transformation from 
u(i) to v{t) = dx/dt. Combining together Exps. (15), (47 1, and (52) we can write this transformation as 
follows 

^=sinh|^+A„signKt)]| (55) 



which is completed by the second proportionality of Exps. (52) relating t to t. 
Summarizing aforesaid we get the following statement. 

Proposition 4. The Levy flights {x{t)} governed by model (pi) with the kinetic coefficients (l7| are equiv- 
alently described by the parameterless stochastic process {u(t)} obeying the stochastic differential equa- 
tion (53 1 with additive white noise. The variable u and the initial variable x are related via Exp. (55) and 



the second proportionality of Exps. (52) 



Propositions [T}[4| are the basic results of the present paper. In what follows we will make use of them 
to demonstrate the fact that model (Isl) does describe the Levy flights of superballistic, quasiballistic, and 
superdiffusive regimes matching < a < 1, a = 1, and 1 < a < 2, respectively. It should be reminded 
that this correspondence was strictly proved only for the superdiffusive Levy flights [3H] and for the other 
regimes it was demonstrated numerically (33] . A more sophisticated analysis of such Levy flights based on 
the constructed representation is worthy of an individual publication. 
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4-2. Asymptotic properties of the core stochastic process 

Let us analyze the probabilistic properties of the peak pattern Pu(A;t) for the core stochastic process u(t) 
in the limit 6^1. To do this, the asymptotic behavior of distribution functions describing the statistics 
of time patterns related to Pu(fct) as —>■ oo will be considered in detail. Dealing with the core stochastic 
process we may set the first boundary U; of the layer L^ equal to zero, u; = 0, without loss of generality. 
The stationary distribution function p^^ (u) of the random variable u is determined by the expression 

p"'(u) == e-" (56) 



after merging the half-spaces {u < 0} and {u > 0}, which stems directly from Exps. (161 and ( [20| ) after 
setting 4'{ri) = 1 and a = 1 in these formulas. 

The probability ©u(t, 0) is determined by the corresponding path integral over all the trajectories {u(t')}o 
meeting the inequality < u(t') < for < t' < t; now the value 9 is noted directly in the list of the 
arguments of the function f&u{l, 9). Therefore the function 

m0):=^^=te-^.exp{-te-«} (57) 

is the probability density that the random variable u attains the maximal value equal to 9, i.e., 

9 = max u(t') (58) 

0<t'<t 

inside the interval (0, t). The latter equality in Exp. (57) stems directly from Propositionlsland formula (56). 
According to Proposition IT] the function C5u (t, 9) admits the interpretation in terms of the peak pattern 
Pu(fct), where the number of peaks or, what is the same, the number of basic units ki := kt = i/(a^r) is 
regarded as a fixed parameter of the random walk classification. In contrast, the specific time moments {ii} 
of the pattern partition including the terminal point t„ = t (Fig. l3| are treated as internal classification 
parameters over which the averaging must be performed. Under these conditions the function 25u(t, 9) is just 
the the probability of the peak pattern Pu(fct) generated by the random variable u and containing exactly 



ki peaks. Therefore the function ^{i,9) introduced via Exp. (57) gives the probability density that the 
maximum attained by the variable u within at least one of the ki peaks is equal to 9. Naturally, the maxima 
attained by this variable inside the other peaks must be less than or equal to 9. Going in a similar way the 
probability density f{9) of the extreme value 9 for one peak is written as 

f(^)-^-ue-. (59) 

Here the value 9 has the meaning 

9= max u(t') (60) 

t'fEa given peak 

and, as before, t^ is the mean duration of individual peaks. Function (57) is plotted in Fig. [6] and its 
asymptotic behavior as — > oo enables us to introduce a single-peak approximation. 

4-2.1. Single-peak approximation 

In the region 9 > 9i, where i > and 

9i = ln(t) (61) 

so that te^^' = 1, the asymptotic behavior of the extreme value probability 5'(0,t) is specified by the 
expression 

Si9,i)^ie-' (62) 



by virtue of (57). Formula (62) can be reproduced in the following way. Let us assume, at first, that the 



maximal value 9 is attained inside a certain peak i. Then inside the other (fct— 1) peaks the random variable 
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Figure 6: Probability density S{6,i) of the extreme value 8 being attained by the random variable u within the peak pattern 
IPuCfci) (curve 1) and its form in the frameworks of the single-peak approximation (curve 2). In plotting these data the time 
t = 50 was used. 



u has to belong to the interval u G {0,9). The probability of the first event is given by the function f{9). 
The probability of the second event can be written as 



9ior 



fet-i 



[i-m] 



fci 



In deriving this formula Exps. (46), ( [56| , and (59) have been used as well as the inequality fct = t/ru ^ 1 
has been taken into account. The latter inequality enables us to confine out consideration to the patten 
configurations for which f{6) ^ 1 and, thus, to ignore the difference between fct and {kt — 1). Therefore, the 
probability of the extreme value 9 being attained inside peak i is 

The choice of peak i is arbitrary. So to find the probability of the extreme value 9 being attained somewhere 
inside the peak pattern Pu(fct) we have to multiply the last expression by the number fct = I/tu of peaks, 
which gives rise to the equality 

mt)-fctfWexp{-fctf(0)} (63) 

coinciding with Exp. (57) by virtue of (59). In the case 9 > 9i the inequality kif{9) ^ 1 holds and the 
exponential multiplier in Exp. (63) can be ignored, leading to Exp. (62). In other words, the extreme value 



9 is attained actually within one peak whereas the effect of the other peaks on its probabilistic properties is 
ignorable. 

Summarizing these arguments we draw the conclusion below. 



Proposition 5 (The single-peak approximation). In the region i 



> 



the probability density ^{9, t) of the 



variable u attaining the extremal value 9 within the peak pattern Pu(fct) exhibits the asymptotic behavior 
that can be represented in the following way. The value 9 is attained exactly within one peak. Variations 
of the random variable u inside the other peaks may be assumed to be small in comparison with 9 and the 
condition u < 9 does not really affect their probabilistic properties. 

Naturally, when a given value of 9 falls outside the region 9 > 9i, i.e., 9 <9i variations of the random variable 
u in many peaks become comparable with 9. Ks a, result, the fact that the condition \x < 9 is imposed on all 
the peaks of Pu(fct) is essential and the single-peak approximation does not hold. In particular, for 9 < 9i 
the probability density ^{9,i) has to deviate substantially from its asymptotics (62) as illustrated in Fig. [6J 
Now let us apply these constructions to the analysis of the particle displacement in the space Mj. . 
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5. Asymptotic properties and scaling law of spatial particle displacement 

The purpose of this section is to demonstrate the fact that for a particle whose random motion is governed 
by model ([3]) the asymptotic behavior of the distribution function 'P{x, t) of its spatial displacement x during 
the time interval t is determined by the asymptotic properties of the peak pattern Pu(fct). By virtue of 
Proposition [5] it can be reformulated as follows. 

Proposition 6. Large fluctuations in the spatial displacement x of the wandering particle during the time 
interval i 3> 1 are related mainly to its motion inside the single peak of the pattern Pu(^t) whose amplitude 
is maximal in comparison with the other peaks and attains extremely large values. 

Proposition [6] is actually the implementation of the single-peak approximation applied to describing the 
particle motion in the space R^- Appealing to the notion of the core stochastic process {u(t)} we also can 
formulate the next statements. 

Proposition 7. The region ^^{t) = {x : \x\ ^ ^{t)} of large fluctuations in the particle spatial displacement 
X directly matches the region > 6*4 of the extreme values 9 attained by the random variable u(t) inside the 
peak pattern Vu{ki)- 

Proposition 8. The asymptotic behavior of the distribution function V{x,t) in the region La;(t) is practi- 
cally specified by the asymptotic behavior of the probability density ^{6, t) inside the region 9 > 9^. 

Proposition 9. The lower boundary x{t) of the region l^x{t) quantifying the characteristic displacement of 
the wandering particle during the time interval t can be evaluated using the single-peak approximation and 
setting 9 equal to the lower boundary of the 0-region wherein the single-peak approximation holds, 9 ^ 9t. 



Below in this section, at first, we will accept Propositions |6]-[9] as hypotheses and using them construct 
the distribution function V{x,t). Then the comparison of the results to be obtained below with the previous 
rigorous results [38l ESI justifies these Propositions. 



Integrating expression (55 1 over the time interval (0, t) we get the relationship between the particle spatial 
displacement x(t) during the given time interval and the time pattern {u(i')}J,^Q of the core stochastic process 

t 

x{t)^^ f sinh I "^ + A„ sign [u(t')] | di' . (64) 



Here the time scale transformation (52) has been taken into account. According to Proposition [6] the time 
integration in Exp. ( 64 ) can be reduced to the integration over the largest peak Fg of the pattern {u(i')}J, Jq 
provided extreme fluctuations of the random variable x{t) are under consideration. In addition. Proposition^ 
claims that the random variable u inside the peak ¥g attains values much larger than unity. The latter holds, 
at least, within a certain neighborhood of the maximum 9 attained by the variable u inside the peak Vg. In 



this case formula ( 64 ) can be rewritten as 

In obtaining this expression we implicitly have assumed without loss of generality that the variable u(t) 
takes positive values inside the peak ¥g. 

The spatial displacement x of the wandering particle and its velocity maximum 

i?=^exp|^| (66) 

attained inside the peak Fg are partly independent variables. Indeed, for example, their ratio 

X 1 f (u(i)-9} , 

- = ^ / exp -^ di 67 
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Figure 7: Example of random trajectories {u(t)} implementing the peak Pg. Initially, t = 0, both of them start from the point 
Uu and terminate when crossing the lower boundary uj = for the first time. This figure also illustrates the technique of 
constructing the distribution function of the maximal value 9 attained by continuous random walks {u(t)} via counting the 
trajectories passing through the layer (9, d + dd) without touching its upper boundary « = 6 + dd. The division of the result 
by the thickness dd of the analyzed layer yields the probability density f(6). 



depends on the details of the pattern u(t) in the vicinity of its maximum 6. Nevertheless, these details seem 
not to be too essential; they determine mainly some cofactors of order unity, see also Ref. |37j . 

To justify the latter statements, first. Figure It] depicts two trajectories u(t) implementing the peak 
P^. The difference in their forms explains the partial independence of the variable x and d. Second, 
Figure Is] demonstrates the statistical properties of such trajectories. The shown patterns were obtained in 
the following way. A collection of random trajectories similar to ones shown in Fig. [7] were generated based 
on Eq. ( [53| with the discretization time step of 0.01. All the trajectories started from the point u„ = 1 and 
terminated when crossing the boundary u; = for the first time. Only the trajectories that passed through 
the layer (6*, 6* + 1) with 6' = 10 without touching the upper boundary + 1 = 11 were taken into account. 
Then for each trajectory the time moment tmax of attaining the corresponding maximum Umax was fixed 
and the trajectory as a whole was shifted along the time axis that the point tmax be located at the time 
origin i = 0. In this way all the trajectories were rearranged that their maxima be located at the same point 
on the time axis. The total number of the trajectories constructed in this way was equal to 10^. Then the 
plane {i, u} was partitioned into cells of 0.1 x 0.1 size and the discretization points of individual trajectories 
fell into each cell were counted. Finally their numbers were renormalized to the obtained maximum. The 
left window in Fig. |8] exhibits the obtained distribution of these values called the distribution pattern of 
u(t) . Actually this pattern visualizes the regular trend in the dynamics of the variable u(t) near the extreme 
point 9 and its scattering around it. As should be expected, the regular trend of u(t) matches the optimal 
trajectory 



'^opt 



{i) = e-\i-u 



of the system motion towards the maximum 



and away from it. The symmetry of this pattern is worthy 
of being noted because only the left branch of the optimal trajectory (68) matching the motion towards the 
maximum (t < tmax) is related to extreme fluctuations in the time dynamics of u(t). The right one is no 
more than a "free" motion of particle under the regular drift. The right window depicts actually the same 
pattern in units of the particle elementary displacement, see Exp. 



(65) 



m, 



6x oc exp < > di 

[ a J 

during the time step dt. It plots this pattern normalized to the particle velocity -d attained at u = 6'. As seen 

19 



distribution pattern of u(t) 



displacement distribution pattern 




0.5 0.0 

time, t 



Figure 8: Spatial patterns visualizing the regular trend in the dynamics of the random variable u(t) in the vicinity of the 
attained maximal value d and the scattering of u(t) around the regular trend. The left window depicts this pattern on the 
plane {t, u}, the right window maps this pattern on the plane {t, 5x}. Here 6x is the elementary displacement of the wandering 
particle along the axis x during the time interval di or, what is actually the same after the correspo ndin g normaliz ation , the 
particle velocity v normalized to the velocity i? corresponding to u = 6. In obtaining these data Eq. ( |53[ l and Exp. | |55| l with 
a = 1 were used; the extreme value was set equal to = 10. The details of constructing the given patterns are described in 
the text. 



in this figure, fluctuations in tlic variable x should be coniparable with its mean value or less than it. In 
constructing the given patterns the value of a = 1 was used. Therefore, in spite of the partial independence 
of the random variables x and d solely the statistical properties of the variable d are responsible for the 
Levy characteristics of the generated random walks. 

For the sake of simplicity in the present analysis we confine our consideration to the regular model of 
the time variations u(t) near the extremum point, i.e., the ansatz 



u(t) 



■^opt 



(t) 



t-tn 



(69) 



will be used in calculating integral ( 65 ) . An approach enabling us to go beyond this approximation will be 



published somewhere else. Substituting ( 69 ) into ( 65 1 we get 



exp 



a 



(70) 



Then using formula (50) and Exp. (62 1 for the extreme value probability density 5^(0, t) we obtain the 
expression 

4ar(l + a) 



Vix,t)^ 



t 



a"r2(a/2) xi+" 



(71) 



giving us the asymptotics of the distribution function V{x,t) of the particle spatial displacement x during 
the time interval t when x 3> x{t) and the expression 



x{t) 



4r(l + a) 
a"r2(a/2) 



-1 1/q 



(72) 
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evaluating the characteristic spatial distance x{t) passed by the particle during the time interval t. Expres- 

nship between th 



sion ( 71 ) has been derived via the relationship between the probability functions 

-1 



dx 

de 



and Exp. (72) was obtained setting 9 = 9^ in Exp. (70 1. 



The rigorous formula for the asymptotic behavior of the function V{x,t) was obtained in Ref. [38j using 
a singular perturbation technique for 1 < a < 2 and then verified numerically also for < a < 1 |39j . 
Following Ref. [55] we rewrite it as 



Whence we see that the rigorous expression and the expression obtained using ansatz ( 69 1 coincide with 
each other within the factor 

'2~ "■ 



n{a)=i-\ r(l + a) e (1,2.05) (74) 

for a G (0,2). The fact that the obtained coefficient ft{a) is really about unity for all the values of the 
parameter a under consideration justifies Propositions [6}|9] 

6. Conclusion and closing remarks 

The work has been devoted to the relationship between the continuous Markovian model for Levy flights 
developed previously |37H40j and their equivalent representation in terms of discrete steps of a wandering 
particle. The present analysis has been confined to the one-dimensional model of continuous random motion 
of a particle with inertia. Its dynamics is studied in terms of random motion on the phase plane {x, v} 
comprising the position x and velocity v — dx/dt of the given particle. Time variations in the particle 
velocity are considered to be governed by a stochastic differential equation whose regular term describes 
"viscous" friction with, maybe, a nonlinear friction coefficient k{v). Its stochastic term containing white 
Gaussian noise, the random Langevin force, allows for the stochastic self-acceleration phenomenon which 
can be of different nature. The stochastic self-acceleration is taken into account via the noise intensity 
g{v) growing with the particle velocity v. Spacial attention is payed to the ideal case where the friction 
coefficient k is constant and the noise intensity g{v) ex v becomes proportional to the particle velocity v 
when the latter exceeds some threshold Va- It is the case where the generated random walks exhibit the 
main properties of Levy flights. Namely, first, the distribution function V{x,t) of the particle displacement 
X during the time interval t possesses the power-law asymptotics. Second, the characteristic length x(t) of 
particle displacement during the time interval t scales with t also according to the power-law. 

The characteristic feature of the considered stochastic process is the fact that such nonlinear dependence 
of the noise intensity on the particle velocity gives rise to a multiscale time pattern {v{t)} and the spatial 
particle displacement is mainly caused by the velocity extreme fluctuations, i.e., large peaks of the given 
pattern [37 . In particular, if we consider the velocity pattern {v{t)} of duration t then the particle displace- 
ment X within the corresponding time interval can be evaluated as x ^ i^tT, where i^t is the amplitude of the 
largest peak available in the given pattern and r is a "microscopic" time scale characterizing the velocity 
correlations. As a result, the statistical properties of the velocity fluctuations cause the Levy time scaling 
of the characteristic length x{t) of particle displacement and endow the distribution function V{x,t) with 
the appropriate power-law asymptotics. 

This feature has made it attractive to represent a trajectory of the wandering particle or, speaking more 
strictly, the time pattern {v{t)} as a sequence of peaks of duration about t and to consider each peak 
as a certain implementation of one discrete step of the particle motion. Unfortunately such an approach 
cannot be constructed directly because the time pattern {v{t)} as a random trajectory is not a smooth 
curve. To overcome this obstacle a complex neighborhood C of the line w = on the phase plane {x, v} 
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has been introduced. It contains two boundaries \v\ = u„ and \v\ — vi, where the choice of the parameters 
Vu, vi meeting the inequahty vi < Uu is determined by the simphcity of mathematical constructions. The 
presence of the two boundaries has enabled us to introduce the notion of random walks outside and inside 
the neighborhood £. The notion of random walks outside C describes a fragment of the particle motion in 
the region \v\ > vi without touching the lower boundary \v\ = vi until the particle gets it for the first time. 
Random walks inside the neighborhood C correspond to a fragment of the particle motion inside the region 
\v\ < Vu without touching the upper boundary |ti| — Vu again until the particle gets it for the first time. For 
the analyzed phenomena the initial particle velocity does not matter provided it is not too large. Therefore 
the initial particle velocity was set equal to v^ without loss of generality. In this way any trajectory of the 
particle motion is represented as a sequence of alternate fragments of random walks inside and outside the 
neighborhood C. A complex unit (basic unit) made of two succeeding fragments of the particle random 
motion inside and outside C may be treated as a continuous implementation of one step of the equivalent 
discrete random walks. The individual duration and the resulting length of these basic units are partly 
correlated random variables. It enables us to regard the constructed representation of random trajectories 
as a certain generalization of continuous time random walks (CTRW). The main difference between the 
CTRW model and the model developed here is the fact that the particle is not assumed to move uniformly 
along the straight line connecting the terminal points of one step. 

For the analyzed model (Is]) it has been demonstrated that the particle motion inside the neighborhood 
C practically contributes only to the duration of the basic units. The particle motion outside C determines 
the spatial displacement as well as contributes to the basic unit duration too. In the given model the kinetic 
coefficients, i.e., the friction coefficient k{v) and the noise intensity g(v) exhibit no singularities at u = 0. 
Therefore the distribution of the basic unit duration has no anomalous properties and, thereby, there is a 
linear relationship 

k — ^t 

between the running time t and the number of basic units fc, naturally, for t ^ t and, so, for fc 3> 1. The 
proportionality coefficient 7 has been obtained as a certain function of the model parameters. As a result, in 
describing statistical properties of such random walks the system may be characterized by the number k of 
basic units imposing no requirements on the duration of the pattern {v{t)} as a whole. It should be pointed 
out that the developed classification of random trajectories holds within rather general assumptions about 
Markovian stochastic processes. So it can be generalized to models where the kinetic coefficients exhibit 
essential singularities in the region of small velocities. However, in this case there is no direct proportionality 
between t and /c, moreover, the integer k has to be treated as a certain random variable; a similar situation 
is met in modeling grain boundary diffusion as a stochastic processes of the subdiffusion type [13 SH] • 

Using the constructed trajectory classification the analyzed stochastic process has been reduced to a 
certain universal stochastic process {u(t)} for which the corresponding governing equation contains no pa- 
rameters. Moreover, the white Gaussian noise enters this equation in the additive manner and the regular 
drift is a piece-wise constant function, namely, the sign(u)-function. It has been called the core stochastic 
process. All the basic parameters characterizing the generated random walks, e.g., the exponent of the Levy 
scaling law, are specified by the model parameters via the corresponding coefficients in the transformation 
from the variable u of the core stochastic process {u(t)} to the original particle velocity v. Namely, it is the 
transformation u — > w as well as the linear transformation of the core process time to the "physical" time, 
t — )■ i. This, in particular, explains us why the main results rigorously obtained for the superdiffusive regime 
of Levy flights hold also for the other possible regimes as demonstrated numerically [3S] . 

In order to elucidate the basic properties of the random walks at hand the developed technique has 
been applied to the core stochastic process to construct the peak pattern Pu(fc) := {u(i)} subjected to 
the condition |u(t)| < 6 for all the moments of time, where 6* ^ 1 is a certain given value. Based on the 
probabilistic properties of the core stochastic process it has enabled us also to introduce the probability 
density S{9,i) that the maximum attained by the variable |u(t)| inside the pattern Pu(fc) is equal to 9. The 
probability density ^{6,i) has been studied in detail in the present paper. In particular, first, it has been 
found out that the region of the asymptotic behavior of ^{0,i) is specified by the inequality ^ ^t, where 
f?t — ln(t). Second, the so-called single-peak approximation has been justified. It states that if the value 
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6 > 9i is attained within a given peak of the pattern Pu(fc) then the variations of the variable u inside the 
other peaks may be assumed to be smaU in comparison with 9 and the condition |u| < 9 does not really affect 
their statistics. Then, as been demonstrated, the asymptotic properties of the generated random walks can 
be formulated in terms of the core stochastic process as follows. 

• On time scales i 3> t large fluctuations in the spatial displacement x of the wandering particle are 
implemented mainly via its motion within the single peak of the pattern Pu(fc) whose amplitude is 
maximal in comparison with the other peaks and attains extremely large values. 

• The region |a;| S> x{t) of large fluctuations in the particle spatial displacement x matches directly the 
region 6* > 6*4 of the extreme values 9 attained by the variable u(t) inside the pattern Pu(fc). 

• The asymptotic behavior of the distribution function V{x, t) in the region |a;| ^ x{t) is practically 
specified by the asymptotic behavior of the probability density 5^(0, i) inside the region 9 > 9i. 

• The lower boundary x{t) of the asymptotic behavior oiV{x^t) determines the characteristic displace- 
ment of the wandering particle during the time interval t. So it can be evaluated using the single-peak 
approximation and setting the value 9 equal 0t, i.e., 9 ^ 9t. It should be reminded that 9t is also the 
lower boundary of the 0-region wherein the single-peak approximation holds. 

In addition it is worthwhile to note that, first, in studying the extreme characteristics of particle motion 
outside the neighborhood C the shape of the velocity pattern {v{t)} can be approximated using the most 
probable trajectory {Uopt(f)}. In particular using these results the asymptotics of the distribution function 
V{x,t) has be constructed and demonstrated to coincide with the rigorous results within a cofactor about 
unity. It opens a gate to modeling such processes in heterogeneous media constructing the most optimal 
trajectories of particle motion in nonuniform environment. 

Second, the first item above can be regarded as a certain implementation of the single-peak approxima- 
tion. It is based on the use of the condition |u(t)| < 6* in constructing the pattern Pu(fc). However, the 
developed classification of random trajectories admits a more sophisticated analysis of anomalous stochastic 
processes. To do this it is necessary to consider a more complex system of restrictions |u(t)|i < 9i with the 
independent external boundaries {9i} for different peaks. In this case it would be possible to speak about 
several peaks of different large amplitude and to go beyond the single-peak approximation. 
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Appendix A. Probabilistic properties of random walks inside and outside the layer C 

It should be noted beforehand that in the present Appendix no approximation of the potential ^{r]) 



introduced by expression (16) will be used, only the general properties (pi) of the kinetic coefficients k{v) 
and g{v) are taken into account. It enables us to make use of the results to be obtained here in further 
generalizations, e.g., to allow for the cutoff effects. 

Appendix A.l. Terminal fragments: the limit case s — >■ and — > 00 

The limit 9 -^ 00 describes the situation when the upper boundary rj = 9 of the analyzed region [0, 9) 3 rj 
is placed rather far away from the origin 7/ = and its effect on the random particle motion is ignorable. In 
this case, from the general point of view, the terminal fragments shown in Fig. [2] are no more than random 
walks starting from at a given point ryo and reaching another point rj in a time t without touching a certain 
boundary r] — C,. Their probabilistic properties are described, in particular, by the probability density 
Q{ri, i|77o, C) of finding the random walker at the point 77 in the time t. The Laplace transform G{'q, s|?7o, C) 
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of this function obeys the fohowing forward Fokker-Planck equation matching the Langevin equation ( 13 1 
(see, e.g., Ref. |50|) 

-dG d$(r;)^ 

drj 



sG 



d_ 
drj 



drj 



8{i] ~ rjo) 



(A.l) 



For random walks inside the layer Cq = [0, C) with the initial point rjQ < ( Eq. (A.l I should be subjected to 
the boundary conditions 



dG d$ ■ 

drj drj 



0, 



GU 



77=0 



r,=C 



0. 



For random walks outside the layer C(^ with rjo > ( the corresponding boundary conditions are 



dG d$ ■ 

drj drj 



0. 



G\ 



0. 



(A.2a) 



(A.2b) 



Since the effect of time in the analyzed phenomena is mainly caused by the properties of the peak pattern 
P(i|fc) and the time scales i ^ 1 are of the primary interest, dealing with the terminal fragments we may 
confine our consideration to the limit s — )■ 0. In this case Eq. (A.l I is reduced to the following 

d IdG d^rj) 
drj drj dij 



-5{Tj-rjQ) 



(A.3) 



and after simple mathematical manipulations using the method of variation of constants for solving difference 
equations we get the desired expression for random walks inside the layer Cq 



C 



G'"(77ho,C) := G(77,sho,C)|r,„<c = £-"*(") 



3"*(''')eH(r?'-77)d??' 



(A.4a) 



and for random walks outside the layer £^ 



Vo 



s-i-O J 

C 



(A.4b) 



where 0h(- ■ •) is the Heaviside step function determined by Exp. (28|. 

The first terminal fragment GJf (77) (Fig. ^ matches the analyzed random walks starting at the point 
rji and reaching the point rj in the time At — t — tn without touching the boundary r/ = ?]„ > rji. So the 
probabilistic weight G™(rj, At) of this fragment is specified by the expression 

g''\fj,At)=girj,At\rjo,C)\vo=m 
and its Laplace transform G"^{rj, s) in the limit s — > takes the form 



G'"(?7,s)| „ = £-"*(") 



Is-J-O 



5"*(V)e(^'„,^)rf^' 



(A.5a) 



by virtue of (A.4a). The second terminal fragment G'^^{rj) (Fig. l2| is also represented by the given random 
walks starting at the point 77„ and reaching the point 77 in the time At = t — tn without touching the 
boundary rji. Thereby the Laplace transform G°^^{rj, s) of its probabilistic weight Q°^^{rj, At) is specified by 
Exp. (A.4b), namely, 



G°"*(77, s)L_,o = e-"*^'') / e"*(''')e(r/ - rj') drj' . 



(A.5b) 
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It should be noted that in both of Exps. (A. 5) the quantity rj can take any arbitrary vahie ry e (0, oo) rather 



than a value from the corresponding interval only. Indeed, if the taken value falls outside this interval the 
relevant expression will give out the probability density equal to zero. 

The latter feature enables us to represent the construction reproducing actually the right-hand side of 



Exp. (26) as 






a*(r,) / ea*('7')rf,^' . (A.6a) 



In deriving this expression the identity Q(r] — r]') + O^rj' — ij) = 1 has been taken into account. In the cause 
under consideration it is assumed that the potential ^{rj) can be approximated by a linear function of rj in 
the region rj ^ 1, namely, ^{rj) ~ C + rj, where C is some constant. Under such conditions the previous 
expression is reduced to 



euiri~m)G°''\v,s)\^^„ + eHiriu-v)G'"'iv,s)\ 



= _g-Q*(»?) 

a 



gQ$(?)„) _ gQ$(j)i) 



(A.6b) 



In particular, these expressions demonstrate that in the limit s — >■ 0, i.e., when t — )■ oo the cumulative 
contribution of the two terminal fragments to the probability density of finding the random walker at 
the point rj is equal to the steady state distribution of the random variable rj within a certain factor of 
proportionality. 

Appendix A. 2. Fragments of the peak pattern V{t \ k) 

Generally the individual fragments of the peak pattern represented in Fig. [3] can be regarded as random 
walks of a particle that initially (i = 0) is located at a point 77 and gets a boundary C, for the first time at 
a moment t. If the initial point 77 is located outside the interval Cq :— [0,0, i.e., rj > C,, then the addition 
condition is imposed on the particle motion; it is not allowed for it to touch or cross the distant boundary 
0^1. The probabilistic properties of these fragments are described by the probability density F{C,,t\ri,0) 
that the particle starting from the point r/ at i = gets the boundary C, at the moment t for the first time 
and, in addition when applicable, during its motion never crosses the distant boundary 9. 



The Laplace transform F{C,, sjyy, 9) of this function determined by Exp. (30) obeys the following backward 



Fokker-Planck equation matching the Langevin equation (131 (see, e.g., Ref. |50j ) 



sF ■ 



d^F d^{r)) dF 



drf 



drj drj 



(A.7) 



For the random walks inside the layer Cq, i.e. when < 77 < C Eq. (A.7 1 should be subjected to the boundary 
conditions 



F\ ,^l. 



dF 
drj 



= 0. 



(A. 



?)=0 



and for the random walks outside the layer Cq, i.e. when (^ < rj < 9 the relevant boundary conditions are 



l,,=C 



= 1. 



F 



= 0. 



(A.8b) 



Since the details of solving Eq. (A.7 1 for the random walks inside and outside the layer >C^ are different we 
will analyze the two cases individually assuming the inequality s ^ 1 to hold beforehand, which matches 
large time scales f S> 1. 
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Appendix A. 2.1. Random walks inside the layer Cq 

Because values of r/; and rju about unity are of the prime interest, see the corresponding discussion in 



Sec. 2.3 here we may consider the thickness of the layer Cq to be also about unity, C '^^ 1- In this case at 
the zero-th approximation in s the solution of Eq. (A. 7) subject to the boundary conditions (A.8a) is equal 
to unity. So at the first approximation in s Eq. (A.7) can be rewritten as 



d'^F d<^{ii) dF 



drf^ 



di] drj 



(A.9) 



Using the method of variation of constants Eq. (A.9 1 under conditions (A.8a) is integrated directly, yielding 
the desired expression 

C r,' 

^'"(C, sh) = 1 - '"^ / */ e"*^"'^ / d?7"e~"*(''"^ . (A.IO) 

-q 

Here the parameter 9 has been omitted from the list of arguments and the superscript 'in' has been added 
to denote the analyzed region explicitly. 

Expression ( |A.10 ) immediately enables us to write the Laplace transform of the individual probabilistic 
weights J^'"(At) of the fragments Fj"^_^i := F'"(At) with Ai = i^+i - t, (Fig- [sj) in the form 



(A.ll) 



F''^(s) = l-s /" d?7' e"*("') /"dry"e-"*("") 

Vi 

within the approximation of large time scales i 3> 1 or, equivalently, for s ^ 1. 

Appendix A. 2. 2. Random, walks outside the layer Cq 

Since the analyzed phenomena are governed by large fluctuations in the particle velocity, the upper 
boundary of the region under consideration 77 S (Ci^) is presumed to be a rather distant point, ^ 1, in 
addition to the assumption s ^ 1. Then let us seek the desired solution of Eq. (A.7 1 subject to the boundary 



conditions (A.8b) in the form 



F(ry)=A,F,(77)+^i^oW, (A.12) 

where As, Aq are some constants and the functions Fs{ri), Fq^tj) are specified via the expression 



Psfiiv) = exp < / ksfl{ri')dri' 



(A.13) 



and, then, to rewrite solution (A.12) in the form 



The boundary conditions (A. 8b) imposed on solution (A.12) enable us to calculate the coefficients As, Aq, 

FsiO) 



^^>~ Fo{9)-Fs{9) ^ '^^> Fo{e)-Fs{e) ^'^^' '^^^^ 



(A.14) 



In deriving Exps. (A.14| the identity Fo{C) — ^s(C) ~ 1 stemming directly from definition (A.13) has been 
taken into account. 



The substitution of (A.13 1 into (A.7) shows the function k{ri) := fcs.o(»7) to obey the Riccati equation 

(A.15) 



drj dr] 
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In the region rj > 1 (when also C ^ 1) the potential ^(jj) introduced via Exp. (161 is approximately a linear 
function of its argument, namely, <&(?y) ~ ry + C, where C is a certain constant, by virtue of the accepted 
assumption (Is]). In this case d^{r])/dr] w 1 and the Riccati equation (A. 15 1 has two solutions 



ksiv) 
ka{Tl) 



s 
a 



(A.16) 

(A.17) 



written in the approximation of leading terms in s ^ 1. Since we can choose any two specific independent 



solutions of the Riccati equation (A. 15 1 let us impose on the desired solutions k — ks.oirj) the requirement 
that in the region r/ > 1 both of them tend to the constant values ( |A.16 ), (A.17), respectively. 

In accordance with the results to be obtained the solution ks{r]) can be treated as a small quantity of 
order s due to the monotonous growth of the potential $(?/) with its argument rj. It enables us to omit the 



quadratic term from the Riccati equation (A. 15) and, then, using the method of variation of constants write 
the desired solution as 



fcsW 



i<i>(j)) 



^-a$(j)') 



drj' . 



(A.18) 



Whence it follows, in particular, that the function Fs{ri) is decreasing one and its asymptotic behavior as 
s — > is specified by the expression 



ri oo 



drj' 



II -a<i>{ri") 



(A.19) 



v' 



For the analyzed phenomena we need to know only the characteristics of small deviations of the function 
Fg (r/) from unity for s <C 1 , which actually has enabled us to construct the function k^ (rj) within the accuracy 
of the leading term in s. 

In constructing the second solution fcg (ry) we may set the left-hand side of equation ( |A.15 ) equal to zero 
and, then, rewrite it as follows 

fc^.^-^. (A.20) 

drj drj 

In the case under consideration the deviation of the function "I>(?y) from the linear dependence can occur 



only on scales 77 ^ 1. So we may confine ourselves to the first iteration in (A.20), yielding us the desired 
approximation 



fco(ry) =a 



d<^{Tj) 



In 



d^{r]) 
drj 



dr/ drj 

Expression (A. 21) and definition (A. 13 1 immediately enable us to specify the fmiction Fo{ri), namely, 

-1 



(A.21) 



Foiv) 



d^{ri) 



dr] 



„a[<E.(,,)-$(C)] 



(A.22) 



In obtaining this formula we have set d^{ri)/dr] = 1 at rj = ( > I. In particular, this formula demonstrates 
us that the function -fo(^/) exhibits the exponential growth as the variable 9 runs to large values. So the 
contribution of the finite magnitude of the parameter s ^ 1 to the function Fq (77) becomes essential only 
when the argument 1] attains such large values that the function Fo{s) even in the leading approximation in 
s, i.e., for s = has already got exponentially large values. This region is of minor interest in the frameworks 
of the present analysis, which has enabled us to construct the solution ^0(77) setting s = 0. 

For the same reasons we may confine ourselves to approximation (A.19) in evaluating the first term on 



the right-hand side of Exp. (A. 14) and in the second term set Fs{r]) = Fs{6) = 1. It is due to the fact that 
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for the most interesting magnitudes of the parameter 6 the effects caused by the finite magnitude of the 

and ( 



parameter s are ignorable. fn this way using formulas (A.f9) and (A. 22 1 we gets the approximation 

V oo 



C n' 



drj 



ri=9 



gQ*(?)) _ ga*(C) 



pa$(6l) 



(A.23) 



Setting here rj — rj^ and C — Vi ^^ immediately get the desired expression for the Laplace transform of the 
individual probabilistic weights J"°"*(Ai) of the fragments ¥°f_^_^ := F°^'(At) with Ai = U+i - U (Fig. |2| 
in the form 






d$(r/) 



drj 



r]=0 



„a*())„) _ pQ*(r;i) 



i$(e) 



(A.24) 



It should be reminded that Exp. (A.24| as well as (A.23) assumes the estimate d^{rj)/dri = 1 to hold for 

Appendix A. 2. 3. Composed element of the peak pattern P(t | k) 

The construction of the peak pattern P(t | fc), see Fig. [s] and Exp. (29), is based on the repetition of the 
same unit, a composed element made of random walks outside the region C followed by random walks inside 
it. The probabilistic weight T{t) of this unit is determined by its Laplace transform equal to the product 



i^(s)=F'"(s).F°"*(s), 



(A.25) 



Whence taking into account formulas (A. 11 1 and (A.24 1 we immediately get the expression 

-gQ*(l)„) _ gQ*(l)i)- 



F{s) = 1 - s / drj'e! 



/ Q*(l)') 



» -a*()) ) 



di' e 



d^{r]) 



drj 



r]=e 



pQ*(e) 



1- 



gQ$(?)„) _ gQ$(j)i) 



rfT?" £-"*("")+ , 



d^{r]) 



drj 



^-a<S>{9) 



(A.26a) 



Using formula ( 20 1 for the stationary distribution function p^* (77) of the random variable 77 and actually 
merging the half-spaces {77 > 0} and {77 < 0} the obtained expression is reduced to the following 



F{s) = 1 - 



^a^ijiu 



„"*('?;) 



dp^*- (77) 



drj 



drj e 



/ -Q*(l)') 



(A.26b) 



It gives us the desired probabilistic weight in the frameworks of the analyzed case, i.e., for large time scales, 
t ^ 1, or equivalently s ^ 1, the distant upper boundary of the region of random walks, 9^1, and the 
accepted assumption about the behavior of the potential ^{rj) in the region 77 '^ 1, namely, $(77) « C + 77, 
where C is a constant. 
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